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ABSTRACT 

Using theoretical formulation of Chung-In Um et. al., the ground 
state energy of two-dimensional liquid 4 He interacting through 
Lennard-Jones and Aziz potential were calculated at various 
densities. Our theoretical result indicates that ground state energy 
of two dimensional liquid 4 He interacting with Aziz potential is 
close to GFMC calculation. 



INTRODUCTION 

In the previous paper, 1 we have 
presented a method of evaluation of radial 
distribution function g(r) of two- 
dimensional liquid 4He interacting through 
the Aziz potential at various densities using 
the theoretical formalism of Chung-In Um 
et. al. 2 In this paper, using the same 
theoretical formalism, we have evaluated the 
ground state energy of two dimensional 
liquid. 4 He interacting through Lennard 
Jones and Aziz potential at various densities. 

For many years, physicist have used 
hypernetted chain (HNC) and Fermi 
hypernetted chain (FHNC) theories to 
investigate the properties of highly 
correlated interacting boson and fermion 

systems such as liquid ^He and normal 

liquid 3 He. The HNC/FHNC equations 



with Jastrow ground-state wave functions 
treat consistently both long and short 
correlations of the system and give 
agreements to some degree with 
experimental results. In spite of qualitative 
successes of the Jastrow ground-state 
function, there exist quantitative differences 
between HNC/FHNC results and 
experimental ones, especially the 
equilibrium ground-state energy, and this 
fact leads one to consider the contributions 
arising from the three-body correlation 
functions. 5 It is well known that the three- 
body correlations do not affect seriously the 
structural properties of ground states, such 
as the radial distribution function and liquid 
structure function, but contribute 

significantly to the ground-state energy. 

3 

Since liquid He system is less dense than 
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liquid He, that is, the equilibrium density of 

the former is 0.0166 while that of the 

latter is 0.02185 A J in a three-dimensional 
system, effects of three-body correlations in 
3 He are rather small compared to those in 

^He. At equilibrium density, they make up 
about half the difference between the 
Jastrow results and experimental energy in a 

, 5 

liquid J He system and lower the Jastrow 
ground-state energy by about 10 % in liquid 
3 He. 

A trial ground-state wave function 
including three-body correlation functions 
in addition to two-body correlation effects 
can be constructed in the form 

y (r,,r 2 ,...r N ) = 

i<i i<j<k 

where f 2 (ry) is the two-body correlation 
function which describes spatial correlations 
and depends only upon the spatial distance 
ry between two particle / 3 (ry, r jk , r Ki ) is a 

three-body correlation function, and I <!>> is 
the anti-symmetric product of non- 
interacting single-particle plane waves 
normalized in the system's area Q. with N 
fermions, which becomes simply unity for 
boson systems. 

In addition to the three-body 
correlation functions in the ground-state 
wave function, one should take into account 
the contributions coming from the 
elementary diagrams as well as effects of the 
HNC/FHNC equations to obtain better 
results. Unfortunately, there does not exist a 
general closed expression for the elementary 



diagrams and therefore one must calculate 
each diagram individually, which is very 
time consuming. During the last two 
decades, two approximations for the 
estimation of the elementary diagrams have 
been developed. One is known as the scaling 

6 

approximation, and the other is the 

7 

interpolating equation approximation. Both 
techniques have shown their efficiency for 
central two-body correlations, providing 
similar results for the total ground-state 
energy to those of the variational Monte 
Carlo simulation. The scaling approximation 
uses the fact that successive contributions of 
the five, six and higher order elementary 
diagrams are approximately proportional to 
and have very similar spatial behaviour to 

6 

the four-body elementary diagrams, so that 
one can represent the total elementary 
diagram contribution by means of the readily 
calculable four-body elementary diagrams. 
The interpolating equation approximation 
determines a parameter such that it gives the 
same value for the classical isothermal 
compressibility obtained from both the 

HNC/FHNC equations and Percus-Yevick 8 
equations using the classical pressure 
derivative or the compressibility integral. 

Although large number of 

4 ^ 
theoretical works regarding He and He 

9-12 

system have been reported but very few 
of them have been used for two-dimensional 
liquid. 

13 

Chung-In Um et. al. have 
developed a theoretical model to study two- 
dimensional quantum fluid. By taking their 
model, one has evaluated the radial 
distribution function g(r) as a function of r 
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for two dimensional liquid He interacting 
through Aziz potential for various densities. 
In section 2, we have given the mathematical 
formulae used in the evaluation. The 
discussion of the result has been given in 
section 3 of this paper. In the last portion, 
we have shown the tables of the results and 
the reference taken. 

MATHEMATICAL FORMULAE USED 
IN THE EVALUATION 

The Hamiltonian for the considered 
system is given by 

H = T + V = XT I+ XV 1J =-^Xvf + Xv(|r 1 -f i |) 

i i<j m i i<j 

(1) 

where V(r) is the interatomic interaction. A 
popular potential used in theoretical 
investigations is that of Lennard-Jones, 
which has two parameters, i.e., the hard- 
core radius and the well depth, and can be 
written as 



V LJ (r) = 4e 



e = 10.22 K, = 2.556 A (2) 
However, Vj^j(r) includes only the 

dipole-dipole interaction and does not take 
into account multiple interactions. A more 
physically realistic potential which accounts 
for the self-consistent field Hartree-Fock 
repulsion and multiple interactions is the 



HFDHE2 potential 
presented as 



of Aziz et. al. 



12 



V u (r) = e' 



A exp -a 



+ c, 



F(r) 



F(r) = 



where 

r exp[-(Dr m /r-l) 2 ] r/r m <D 
1 r/r m <D 
The values of the constants are 
A = 0.54485 x 10 e 



(4) 



e*/k B = 10.* 



a = 13.353384, 
C 8 = 0.4253785, 
D =1.241314, 



10 



l m 



= 1.3732412 
= 0.178100 
= 2.9673 A 



f 0 ^ 


12 




6 


X 













It is well known from the Green's 
function Monte Carlo simulations and other 
variational calculations that the Aziz 
potential gives closer results to experiments 
than Vlj in three dimensions. One uses both 

potentials and compare the results with those 
in two dimensions. 

The n-particle distribution function is 
represented as 

p n =(?„...?„) = v n N(N-l)...(N-l) 
J\|/ + (f 1 ,...f n )\|/(f 1 ,...f n )dr n+1 ...df N 



(4) 



where r t denotes the spatial coordinates of 
the i-th particle and v is the spin degeneracy 
of the system (2 for this system). In the 
homogeneous system, as N — > <*> and £1— >°°, 
the single-particle distribution function 
reduces to the density of the system so that 
the radial distribution function can be 
expressed interms of the density and two- 
particle distribution function P2O") as 



g(r)=p2P 2 (r) = 



(3) 



N(N-1) 



Jv|/ + (f 1 ,...f n )\|/(f 1 ....f n )df 3 ...df N 



(5) 
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One uses the FHNC approximation 
to sum the diagrams arising from the cluster 
properties of g(r). In this scheme, g(r) can be 
decomposed as 



g(r y ) = gdd(ri J ) + 2g de (r 1J ) + g ee (r 1J ) 



(6) 



where dd, de and ee represent terms in 
which both i and j are not exchanged, only j 
is exchanged, and both i and j are 
exchanged, respectively. The components of 
g(r-) are given by 

gdd(rg) = exp[u 2 (ry) + N dd (r y ) + E dd (r y )] 
gde(r y ) = gdd(r y )[N de (r y ) + E de (r y )] ^ 

&efo) =g44-i%) /V+ M%) +E 4 r 9 ) + K( r 9 ) +E M[ 

geeh) = g dd (r y )L(r y )/v 
where 

L(r y ) = -l(k F r y ) + v(N ee (r y ) + E ee (r y )) 

(8) 

l(x) = 2j 1 (x)/x (9) 

k F is the Fermi momentum of the system, 
and ji(x) is the Bessel's function of the first 
kind of order 1. N mm represent sums of the 
nodal diagrams and E nm sums of the 
elementary diagrams. The equation g ee (rij) 
denotes the terms in which both i and j are 
exchanged in an incomplete exchange loop, 
and N^ry) and E ee (ry) are sums of the nodal 
diagrams and the elementary diagrams in 
which i and j belong to the same permutation 
loop, respectively. 

With the FHNC equation and 
expression for the g(r), one now calculates 
the energy per particle of the system 



E var= i (v|%) = <y) , (T) 

N N N(yly) N N 



(10) 



The potential energy per particle is 
expressed by taking the thermodynamic 
limit, i.e., the area of system £2 and the 
number of particles N are infinite, keeping 
the density constant as 



V 



N 



(ii) 



The expression for the or the kinetic 
energy part is more complicated. There are 
three well-known expression for the helium 

system i.e., Jackson-Freeberg, ^ Clark- 

11 1 o 

Westhous 11 and Bethe-Pandharipande. 
Among expressions for the kinetic energy 
per particle. The Jacson-Feenberg form is 
popular and easily tractable since its three- 
particle distribution part is not more 
significant than others. One uses here the 
Jackson-Feenberg and the Bethe- 
Pandharipande forms so that we make use of 
the scaling approximation for the elementary 
diagrams. From the identity 



J \|/ + V 2 Ydf = ^ J [y + (V>)(V V) + -2(VY + ).(VY)dr] 

(12) 



the Jackson-Feenberg energy form is 
obtained as 

~ %2 , 1 V(<i>|F 2 V 2 V 2 +i-FV 2 F--(V.F) 2 |<D) 
N 2mN Y 1 ' 2 ' 2 V 1 n 1 

<V> 
N 

(13) 



2 m N (\|/ I \|/ ) 



^-X F 2 V 2 |4>| 2 dr,...dr N + 



= T F =W B +W,+U^ 
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The first term of Eq. (13), where W and U are the bosonic energies per 
*2i,2 , A „ particle, and Wp and Up are from the last 

term of Eq. (19) 



T F ft k F /4m^ is the non i n t er acting 
fermion energy with the Fermi momentum 
k F - -y/47lp / v ^ an( j ^ second, Wg is the 
energy of the boson system. 



w. 



n =^Jd 2 rg(r) 



V(r)- 



* 2 V 2 / 2 (r)-(V/ 2 (r)) 2 



(2m) 



The forms of and are 



V 2 l(k F r) + Vl(k F r) 2 -vd dd (r)E cc (r)V 2 l(k F r)} 

(15) 



U 4 = 



8m 



J d 2 r l2 d 2 r l3 (g dd (r l2 ) - l)Vl(k F r 12 ).(g dd (r 13 ) - 1) 

xVl(k F r 13 )g ee (r 23 ) 

U x =^ 3 -Jd 2 r 12 d 2 r 13 g3(r 12 ,r 2 3,r 1 3)V 2 u 3 (r 12 ,r 2 3,r 



137 

(17) 



The Pandharipande-Bethe energy 
form is obtained by expanding the kinetic 
energy per particle as 

ft 2 



2mN ; 

which yields 

(T) PB * 



l — y/4)|FV 2 [F|4>)l = ^^ 
\Hf)i 1 L 1 'J N 



, N Y (* |F[V 2 + Vf + 2(V i F)| <D) 

N 2mN (yly) j L v . yj / 

(19) 

After some calculations, we obtained Epg/N 



as 
E 



PB 



N 



= T F +W + W r + U + U F 



W : 



u = 



:jd 2 r g (r 



V(r)- 



*l)v 2 / 2 (r) 2 

m 



/z(r) 



(21) 



ft p 



2m 



Jd 2 r 12 d 2 r 13 g 3 (r 12 r 23 ,r 13 ) 



v./jM-v./A,) 



(22) 



(14) W F =-^Pjd 2 rg K (r) 



2nr 



V/ 2 (r).Vl(k F r) 



(23) 



u F = - 



^Jd 2 r 12 V/ 2 (r 12 ).Vl(k F r 13 ) 



2m 



x{gee(r, 3 )g dy (r 1 2)(g 

dy( r 32 )-l) + gdd(ri 3 )gee(r23)gee(r 3 l)} 

(24) 

For numerical calculations, one 
adopts the McMillan type function as a trial 
(16) two-body correlation functions, which is 
widely used in variational and Monte Carlo 
calculations for liquid and solid helium 
systems. 

u 2 (r) = - 



(25) 



where b is a variational parameter chosen to 
give an optimized radial distribution 
function at a given density. To obtain more 
precise results, we include three-body 



(18) correlation effects in the following forms: 
u 3 (r iJ ,r jk ,r ki ) = ^Ti(rjTi(r jk )(f ij .f jk ) 



r|(r) = 7Vexp 



eye 



r-r, 



V ro i J 



(26) 



(27) 



where eye denotes the cyclic permutation 
among the three particle coordinates i, j and 



(20) k, and 1J is a unit vector along the line 
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connecting particle i and j. Since it was first 
used in a variational Monte Carlo 

13 

calculation, this form for the three-body 
correlations is generally used in variational 

HNC/FHN methods.' 4 



DISCUSSION OF RESULTS 

In this paper, we have evaluated the 
ground state energy of two dimensional 

liquid ^He interacting through Lennard- 
Jones and Aziz potential at various densities. 
We have used the theoretical formalism of 

2 

Chung-In Um et. al. in our evaluation. 
Chung-In Um^ calculated ground state 
energy of two dimensional liquid ^He with a 



Jastrow two body correlation function in the 
trial wavefunction and within the HNC 
approximation. His other calculation for 
energy per particle is by taking two body 
correlation and contribution of four order 
elementary diagrams. We have taken three 
body as well as tow-body correlations 

15 

within the scaling approximation. Our 
theoretical calculation indicates that the two 

dimensional liquid ^He system interesting 
through Aziz potential in more close to the 

14 

GFMC calculation. Comparing with two 
potential, one conclude that Aziz potential is 
more suitable to describe, the ground state of 
the two dimensional helium system, similar 



to bulk case. 
Table T, 



16 



Evaluated results of ground state energy of two-dimensional liquid He interacting through the 
Lennard- Jones potential at various density 



P(A~ 2 ) 


Ground state energy E(k)/N 


Our's 


Chang's results 


GFMC Calculations 14 


0.020 


- 0.520 


- 0.628 


- 0.504 


0.025 


- 0.546 


- 0.637 


- 0.546 


0.030 


- 0.624 


- 0.658 


-0.589 


0.035 


-0.658 


-0.586 


- 0.627 


0.040 


- 0.748 


- 0.552 


- 0.553 


0.045 


-0.786 


- 0.505 


- 0.524 


0.050 


-0.842 


- 0.486 


- 0.508 


0.055 


- 0.725 


- 0.453 


- 0.493 


0.060 


-0.556 


- 0.446 


- 0.486 
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Table T 2 

4 

Evaluated results of ground state energy [E(k)/N] of two-dimensional liquid He interacting 
through Aziz potential at various densities 



P(A~ 2 ) 


Ground state energy E(k)/N 


Our's 


Chang's results 


GFMC Calculations 14 


0.020 


- 0.487 


- 0.586 


- 0.523 


0.025 


- 0.509 


- 0.603 


-0.588 


0.030 


-0.586 


- 0.658 


- 0.602 


0.035 


- 0.627 


- 0.709 


- 0.657 


0.040 


-0.723 


- 0.682 


- 0.697 


0.045 


- 0.706 


- 0.596 


- 0.602 


0.050 


-0.658 


- 0.542 


- 0.524 


0.055 


-0.532 


- 0.506 


- 0.486 


0.060 


- 0.456 


- 0.487 


- 0.442 
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